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ABSTRACT 

Adaptive Optics observations have dramatically improved the quality and versatility of high angular resolution 
measurements of the center of our Galaxy. In this paper, we quantify the quality of our Adaptive Optics 
observations and report on the astrometric precision for the young stellar population that appears to reside in a 
stellar disk structure in the central parsec. We show that with our improved astrometry and a 16 year baseline, 
including 10 years of speckle and 6 years of laser guide star AO imaging, we reliably detect accelerations in the 
plane of the sky as small as 70 /^as yr~ 2 (~2.5 km s _1 yr _1 ) and out to a projected radius from the supermassive 
black hole of 1'.'5 (~0.06 pc). With an increase in sensitivity to accelerations by a factor of ~6 over our previous 
efforts, we are able to directly probe the kinematic structure of the young stellar disk, which appears to have an 
inner radius of 0'.'8. We find that candidate disk members are on eccentric orbits, with a mean eccentricity of 
< e > = 0.30 ± 0.07. Such eccentricities cannot be explained by the relaxation of a circular disk with a normal 
initial mass function, which suggests the existence of a top-heavy IMF or formation in an initially eccentric disk. 
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1. INTRODUCTION 

The proximity of the Galactic center (GC) makes the region uniquely suited for the study of a supermassive black 
hole and its effects on the surrounding nuclear stellar cluster. Advancements in technology, in particular adaptive 
optics (AO) systems on large (8-10 m), ground-based telescopes, have afforded us a view of this dense region at 
such high angular resolution that individual stellar orbits of short-period stars can be precisely measured. The 
orbit of the central arcsecond star, SO-2, for example, has provided the best evidence to date that supermassive 
black holes (SMBH) exist at the centers of normal galaxies, and that the black hole at the Galactic center, Sgr 
A*, has a mass of Mbh = 4.1 x 10 6 M Q and distance of R = 8.0 kpc, respectively!^ 

With high precision astrometric and spectroscopic measurements, much can also be understood about the 
nuclear star cluster in the vicinity of the SMBH. While the GC is populated predominantly by old, late-type 
stars (~1 Gyr), spectroscopic observations of the region have unveiled ^200 massive, young (~ 6 Myr) stars 
in the central parsecP^ whose origins are a puzzle due to the strong tidal forces from the black hole. Given 
the youth of this population and the relaxation timescale in the GC (~1 Gyr) P the stars' dynamical properties 
contain a signature of their formation history. A particularly prominent feature that has been observed is a 
stellar disk containing a significant fraction (possibly up to 50%) of the young stars orbiting the black hole in 
a clockwise senseP^HED The orientation of the disk orbital plane can be described by its inclination and angle 
to the ascending node, and is roughly i ~ 130° and VL ~ 100°, with a la dispersion of ~lO°Pll2]ll2] The disk's 
surface density profile has been observed to fall off like i? -2 , which has been the main argument for in situ 
formation and against cluster 

While there is consensus in the literature regarding the existence of the clockwise disk and its surface density 
profile, many of its properties have yet to be constrained, in part because interpretations of kinematic studies 
rely on the ability to assign disk membership. The difficulty in determining which stars belong to the disk 
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arises from the lack of knowledge about a star's line-of-sight distance relative to Sgr A*. This can be overcome, 
however, with detections of accelerations in the plane of the skyP^ While some stars in the central arcsecond 
(r < 0'.'5) have been seen to orbit Sgr A* completelyP^l the detection of accelerations becomes more difficult at 
larger radii where the clockwise disk stars are located (r > Of! 8) and requires both high precision astrometry and 
patience. In this proceeding, we report on the quality of our AO observations (|J3J), our astrometric precision, 
including detections of accelerations in the plane of the sky for the young stars ($3]), and finally the eccentricity 
distribution of the accelerating disk stars (Q. 

2. HIGH ANGULAR RESOLUTION IMAGING 

The Galactic center region was imaged at high angular resolution with the Keck telescopes from 1995-2011. In 
the first decade, the images were made with the near-infrared speckle imaging camera, which provided diffraction 
limited images (8 — Of! 05) with Strehl ratios of ^0.05. In each final image constructed from the speckle data 
(< N exp > ~ 1600, t exp ~ 0.14 s, coadd =1), we detect ~160 stars on average, down to K ~ 16 in a 5"2x3!2 
field of view. 

Starting in 2004, images were obtained with the Keck II laser guide star adaptive optics system in conjunction 
with the facility near infrared camera, NIRC2 (< N exp > ~ 100, t exp ~ 2.8 s, coadd = 10). The tip/tilt star is 
located at a distance of r = 19" from the GC. Despite this relatively large distance, as compared with speckle 
imaging these data provide nearly an order of magnitude improvement in Strehl ratio (~0.28), a factor of 4 
increase in areal coverage (field of view of 10"xl0") and a considerable increase in the average number of 
detected stars (^2000) down to K ~ 20. Details of these observations can be found in Ghez et al.J^Lu et al.p^ 
and Yelda et al. (in preparation). With these data, we measure the proper motions and accelerations in the 
plane of the sky for a sample of ~70 young stars, which have known line of sight velocities (t^plUMMI an( j which 
are located at projected radii r > 0'.'8 - 7". 

3. HIGH PRECISION ASTROMETRY 

There are several sources of statistical uncertainty in the positional measurements that must be accounted for. 
While speckle images have typical centroiding uncertainties of a pos ~ 0.85 mas, the high angular resolution AO 
images allow for nearly an order of magnitude improvement for the brightest stars with a pos ~ 0.1 mas (red 
squares in Figure [lj. This improvement is a result of the higher signal to noise ratio in the AO images. A second 
source of uncertainty comes from the transformation process that aligns the stars' positions across all epochs 
(a a in)- These alignment errors are a function of time from the reference epoch and of the number of reference 
stars used in the transformation to the so-called "maser-frame" in which Sgr A* is at rest (see Yelda et alP^ for 
details of this process). As seen in Figure [U a a i n is minimized near the reference epoch, 2006 June, and is larger 
for the speckle epochs (a a i n > 0.5 mas), which have on average ~6x fewer reference stars than are available in 
AO epochs. Finally, there is also a source of statistical uncertainty in the positional measurements that is not 
captured in the other terms that arises from inaccuracies in the estimates of the PSF wings of neighboring sources. 
This additive error term (a a dd) is estimated as the amount of error required to minimize the residuals between 
the observed and expected velocity y 2 distribution for a given number of degrees of freedom (see Clarkson et 
alP^ and Yelda et al. (in preparation) for details). The size of this error term is negligible for the speckle data 
(&add = 0.18 mas), but becomes important for the AO images (<J a dd = 0.10 mas). 

Positions as a function of time are fit for proper motions and accelerations. Whether a star exhibits significant 
accelerated motion in the plane of the sky depends on several factors, including its distance from the supermassive 
black hole, the time baseline over which it is detected, and the precision with which its positions are measured. In 
order to determine whether a star's motion is best described by a velocity or acceleration model, we first compare 
the goodness of fit of the two models using the F-test, in which we estimate the quantity F — X pm /Xaccj which 
follows an F-distribution under the null hypothesis.^ Here Xpm an d Xacc are the reduced x 2 values for the 
proper motion and acceleration fits, respectively. We require that F is greater than the critical F value at the 
4ct significance level and that a star is detected in more than 30 out of the 45 possible epochs, which removes 
the sources that show unphysical accelerations arising from effects such as source confusion or edge effects in 
speckle images. The resulting six stars have reliable radial acceleration measurements at >5 sigma. For all other 
sources, the proper motion fit is used. 
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Figure 1. Alignment [filled black points) and centroiding [unfilled red squares) uncertainties as a function of epoch. The 
median uncertainty of the young stars is reported for each epoch. Alignment errors are minimized near the reference 
epoch, 2006 June, and increase with time away from this epoch. All epochs with a a i„ > 0.5 mas are from speckle imaging, 
where the higher uncertainties are a result of very few reference stars as compared to AO data. The additive errors for 
speckle and AO are shown as dashed lines. 



The position, proper motion, and acceleration errors from the fitting procedure as a function of K magnitude 
(A = 2.2 fxm) are shown in Figure [2] for our sample of young stars. The median errors in positions and proper 
motions are 0.05 mas and 0.03 mas yr _1 , respectively. Stars detected in relatively more epochs show smaller 
errors in positions and velocities, as seen in the top panels of Figure O Acceleration errors are plotted for the six 
stars that have reliable acceleration measurements, as well as 11 stars with 3cr upper limits below their theoretical 
maximum acceleration (see below). The average acceleration uncertainty among these stars is 10 //as yr~ 2 (~0.4 
km s _1 yr -1 ), which is a factor of six improvement over our earlier efforts in Lu et alP^ 

Figure [3] shows the radial acceleration measurements as a function of projected radius for stars in the central 
10" field for which we have acceleration information. The six significant accelerations are shown with la error 
bars. This increases the number of acceleration measurements beyond 1" over our previous work in Lu et alP^ 
by a factor of six, or equivalently, an additional five stars, three of which are reported by Gillessen et alP The 
most distant star from the SMBH for which an acceleration measurement is made is IRS 16SW, located at R = 
l'/5 (~0.06 pc), which is well outside the inner edge of the stellar disk. An additional 11 stars have 3er upper 
limits below the theoretical maximum acceleration given by each star's projected radius, thereby requiring a 
non-zero line-of-sight distance. 

4. ORBITS OF THE GALACTIC CENTER YOUNG STARS 

With six kinematic variables measured (xq, yo, v x , v y , v z , a/{), the orbital elements can be estimated if the 
black hole properties are known. We use updated SMBH properties, which are described in detail in Yelda et 
al. (in preparation). We estimate uncertainties in the orbital period (P), eccentricity (e), inclination (i), angle 
to the ascending node (f2), longitude of periapse (uj), and time of periapase passage (T ) by carrying out a 
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Figure 2. Position (top left), proper motion (top right), and radial acceleration errors (bottom) as a function of K magnitude 
for our sample of young stars beyond a projected radius of 0'.'8. Errors are estimated from either the proper motion or 
acceleration fit to each star's individual positions over time. Stars with N > 30 epochs (filled red points) have smaller 
errors in position and proper motion than stars in fewer epochs (open circles). The radial acceleration errors are shown for 
stars in N > 30 epochs, where we separately plot those six stars passing the F test for accelerations (open red squares) and 
those stars for which no detectable accelerations are found (filled black points), requiring a non-zero line-of-sight distance. 
The average acceleration uncertainty among these stars is 10 pas yr" ~ 2 (-0.4 km s" 1 yr" 1 ). 



Monte Carlo simulation in which we create 10 5 artificial data sets. In each data set, we sample the six kinematic 
measurements, as well as the gravitational potential parameters (the black hole's mass, distance, position, and 
an extended mass component). For stars without or, we assume a uniform acceleration prior, with the minimum 
acceleration set by assuming the stars are bound and the maximum acceleration set by the stars' projected 
radii. Radial velocities, v z , are measured with the Keck AO-assisted integral field spectrograph, OSIRISp^or 
taken from the literature if neededP^ The kinematic variables are transformed to six orbital elements using 
the analytic orbit equations presented in Lu et alP^l 




0.04 0.06 0.08 0.10 0.12 
Projected Radius (pc) 

Figure 3. Accelerations along the radial coordinate as a function of projected radius. The theoretical maximum acceleration 
(|a|zo) for the black hole mass of 4.6xl0 6 Mq (see Yelda et al., in preparation) is shown as the dashed curve, with the la 
upper and lower boundaries shown as dotted curves. For our sample of stars at projected radii r > 0'.'8, we reliably detect 
six significant accelerations out to R=l'.'5 (0.06 pc), shown with la error bars. These accelerations yield line-of-sight 
distances to these sources, allowing for precise determination of the stars' orbits. Stars with 3<r acceleration upper limits 
below the theoretical maximum acceleration (+la) are shown as downward pointing arrows and have strong constraints 
on their line-of-sight distances (|z| > 0). 



For the six accelerating stars, their line-of-sight distances are determined to within ^0.01 pc (la) and therefore 
result in precisely determined orbital parameter estimates. Each of these sources has an inclination, i, and angle 
to the ascending node, ft, that is consistent with the orientation of the clockwise diskJ^E^ES] making all six stars 
candidate disk members (Yelda et al., in preparation). Figure 3] shows the combined eccentricity distribution 
for these sources, where only disk solutions (i.e., those solutions that are within 15° from the orientation of the 
disk) are included. The orbits for these stars have a mean eccentricity of < e > = 0.30 ± 0.07. 

5. CONCLUSIONS 

Adaptive optics observations have provided a substantial leap forward in the quality of Galactic center measure- 
ments. As compared to speckle imaging observations, nearly an order of magnitude improvement is seen in Strehl 
ratio, number of stars detected, and astrometric precision. Such high astrometric precision of stars in the GC 
has led to detections of accelerations in the plane of the sky, which thereby give a measurement of line-of-sight 
distances from the black hole. We present the accelerations and eccentricity distribution for six young stars, 
which are all candidate members of the clockwise stellar disk in the central parsec. The mean eccentricity of the 
disk stars is < e > = 0.30 ± 0.07. Given the age of these stars, self-relaxation in an initially circular disk with 
a normal IMF cannot excite the orbital eccentricities to such high values. Alexander et alP^ showed that such 
eccentricities can be reached only if the IMF were extremely top-heavy. Although some authors have claimed 
evidence for a top-heavy IMFpEl the matter is still under debate (Lu et al., in preparation). Alternatively, the 
stars could have formed in an originally eccentric disk. 

Astrometric measurements of the Galactic center arc currently limited by knowledge of the point spread 
function, which has been assumed to be constant across the image. However, both instrumental and anisoplanatic 
effects lead to a spatially- and temporally-variable PSF. Our group is currently developing a methodology to 
accurately model both the static, instrumental PSF (see Fitzgerald et al., this volume) and the time- variable PSF 



5 



4 - 




o 

CL 



1 - 



ft.O 0.2 0.4 0.6 0.8 1.0 

Eccentricity 

Figure 4. Eccentricity distribution from 10 Monte Carlo trials for the six accelerating stars. Solutions that are within 
15° of the clockwise disk orientation are included. All six stars have orbital planes that are consistent with the clockwise 
disk. These candidate disk stars have eccentric orbits, with an average eccentricity of 0.30 ± 0.07. 



using atmospheric profiler data taken simultaneously with our science observations. Accurate PSF modeling will 
allow for an improved optical distortion model and a more stable astrometric reference frame, which is critical 
for detecting the effects of general relativity and extended mass on the orbits of short-period stars at the GC.^ 
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